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ABSTRACT 

We apply magnetohydrodynamic (MHD) modeling to the radio galaxy Hercules A for investigating the jet- 
driven shock, jet/lobe transition, wiggling, and magnetic field distribution associated with this source. The 
model consists of magnetic tower jets in a galaxy cluster environment, which has been discussed in a series of 
our papers. The profile of underlying ambient gas plays an important role in jet-lobe morphology. The balance 
between the magnetic pressure generated by axial current and the ambient gas pressure can determine the lobe 
radius. The jet body is confined jointly by the external pressure and gravity inside the cluster core radius 
R c , while outside R c it expands radially to form fat lobes in a steeply decreasing ambient thermal pressure 
gradient. The current-carrying jets are responsible for generating a strong, tightly wound helical magnetic 
field. This magnetic configuration will be unstable against the current-driven kink mode and it visibly grows 
beyond R c where a separation between the jet forward and return currents occurs. The reversed pinch profile 
of global magnetic field associated with the jet and lobes produces projected 23-vector distributions aligned 
with the jet flow and the lobe edge. AGN-driven shock powered by the expanding magnetic tower jet surrounds 
the jet/lobe structure and heats the ambient ICM. The lobes expand subsonically; no obvious hot spots are 
produced at the heads of lobes. Several key features in our MHD modeling may be qualitatively supported by 
the observations of Hercules A. 

Subject headings: galaxies:individual: Hercules A — galaxies: active — galaxies: jets — methods: numerical 
— MHD 



1. INTRODUCTION 

In the present paper, we continue our discussion of 
the dynamics of extragalactic jets in cluster environments 
within the framework of the "magnetic tow er" model 
dLvnden-Bell & Boily|[T994t lLvnden-Beli1 IT996). which we 
have analyzed in a series of papers. Magnetohydrodynamic 
(MHD) mechanisms are frequently invoked to model the 
launching, acceleration and collimation of a strophysical jets 
(see, e.g., Ferrari 1998: iMeier et al. 1 1200 U and refer ences 
therein). An underlying large-scale (coronal) poloidal field 
for producing the magnetically driven jets is almost univer- 
sally assumed in many theoretical/numerical models. How- 
ever, the origin and existence of such a galactic magnetic field 
are still poorly understood. 

In contrast with the large- s cale field models, L ynden-Bell 
dLvnden-Bell & Boilv I [l994t iLvnden-Belil 1 19961) examined 
the expansion of the local force-free magnetic loops an- 
chored to the star and the accretion disk by using the semi- 
analytic approach. Twisted magnetic fluxes due to the disk 
rotation make the magnetic loops unstable and splay out 
at a semi-angle 60° from the rotational axis of the disk. 
Global magnetostatic solutions of magnetic towers with ex- 
ternal thermal pressure were also computed by iLi et al. I 
using the Grad-Shafranov equation in axisymme- 
try (see also. | Lovelace et al. Il2002|: Lovelace & Romanoval 
120031: lUzdenskv & MacFadvenll2006l) . Full MHD numer- 
ical simulations of magnetic towe rs have been performed 
in two-dimension ( a xisymmetric) dRomanoya et al.l 119981: 
iTurner et al. I 119991: lUstvugova et al. I 120001: iKudoh et al. I 
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2002t iKato et al. 1 12004) and three-dimension (Kato et al. 



2004). Magneti c towers are also observed in labora tory ex- 



periments dHsu & Bellan ll2002tlLebedev et al. 1120051) . 

The first in our series, ILi et al. I (120061 hereafter Paper I), 
described the basic assumptions and approaches in the nu- 
merical modeling of magnetic tower jets. The evolution of 
the tower jets in a constant density /pressure background was 
exami ned there. The second in our series, Nakamur a et al. I 
(2006, hereafter Paper II) investigated the global structure of 
magnetic tower jets in a gravitationally-stratified atmosphere, 
in terms of the MHD wave structure, radial force equilibrium, 
and collimation. A large current flowing parallel to the jet 
bulk flow plays an essential role in the determining the lobe 
radius, as does the background pre ssure profile. The third in 
our series, Nak amura et al. I d2007i hereafter Paper III) inves- 
tigated the stability properties of magnetic tower jets on clus- 
ter scales. Current-driven instabilities are responsible for the 
non-axisymmetric structures; the external "kink" (m = 1) 
mode grows outside R c , while the internal "double helix" 
(m = 2) mode grows predominantly inside R c . This is the 
forth in our series of papers to examine the nonlinear mag- 
netic tower jets. Here we examine our model's applicability 
to an individual source, Hercules (Here) A (3C 348), by per- 
forming three-dimensional MHD simulations. 

Here A is one of the most powerful double-lobed radio 
sources in the sky. Its total energy content is approximately 
3 x 10 60 erg and it is identified with the central comp act 
diffuse (cD) galaxy of a cluster ( S adun & Haves 1 1 19931) at 
a low redshift of z = 0.154 (Sieb ert et al. I Il999l). The 
X-ray luminosity of 4.8 x 10 44 erg s" 1 dGizani & Leahy I 
2004) is thought to be due to bremsstrahlung of the very 
high temperature intracluster medium (ICM). Here A has a 
peculiar radio structure; the morphology looks like a classi- 
cal double-lobed Fanaro ff-Riley type II (FR II) radio source 
dFanaroff & Riley 111 9741) . but it has no compact hotspots, in- 
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stead showing an unusual jet-dominated morphology as in an 
Fanaroff- Riley type I (FR I) radio sou rce, which was first re- 
vealed by Dreher & Feigelson] d!984f) . Furthermore, the two 
jets in Here A are quite different in appearance. The east- 
ern jet appears to have continuous twisting flows, while the 
western jet leads to a unique sequence of "rings". The jets 
in Here A are well-collimated initially but they flare suddenly 
at a certain distance from the nucleus. An explanation for 
the presence of such different structures in the same source 
is still unresolved in the literature. Here A is not a typical 
FR I source, either; its jets are well-collimated and has knot 
components. Therefore, Here A might be classified an inter- 
mediate case: FR I/IT. 

Various models have been proposed for explaining the 
formation of the jets and rings in Here A, including both 



kinematic dMason et al. I [1988; Mor rison & Sadun 
Sadun & Morrison 1 12002|) and dynamic dMeier et al. I 



1996 



1991 



Saxton et al. 1 12002[) points of view. By performing ax- 
isymmetric, two-d imensional hydrodynamic simulations, 
iMeier et al. I dl991l) proposed a model in which an over-dense 
and over-pressured jet undergoes a sudden expansion when it 
becomes highly over-pressured compared to the background 
gas at a certain distance from the core. This causes a subse- 
quent conical expansion, giving rise to the "ice cream cone" 
shape of the lobe. However, at present, there is nothing in the 
X-ray observations to suggest that the jet becom es thermally 
over-p ressured as required in this model. In ISaxton et al. I 
(2002), the formation of ring structure in the western jet is in- 
terpreted as the result of ribbon-like, annular shocks propagat- 
ing through the jet cocoon. Also, these authors conclude that 
the absence of hot spots in Here A results from the dynamics 
of turbulent and entraining ba c kflow s in the radio lobe. 

Recently, Giz ani & Leahy I ([2003) have revealed detailed 
structures of total intensity, spectral index, polarization, and 
projected magnetic field in Here A by using multifrequency 
VLA imaging. The jet disruption in the eastern lobe is clearly 
visible, implying flow instabilities. The projected magnetic 
field, after correcting for Faraday rotation, closely follows 
the edges of the lobes, the jets, and the rings; the two lobes 
have generally similar field patterns. This observational pic- 
ture may disagr ee with the t urbule nt structure inside the lobes 
suggested bv lSaxton et al. I (120021) . 

Magnetic field distributions associated with the extragalac- 
tic jets and their comparisons with the observations have been 
examined by using Stokes parameters, Faraday rotation mea- 
sure (RM), and the projected magn e tic fie l d (-B-vector) (e.g., 
Laingl [19811; IClarke et al. 1 [19891 fl99l lHardee & Rosenl 
1999L 120021: lLaing et al.1l2006h . In general, helical magnetic 



fields produce asymmetric transverse brightness and polariza- 
tion profiles; they are symmetrical only if the field is purel y 
toroidal or the jet i s at 90° to the line of sight (lLaing 111981b . 
Laing et al J d2006l) conclude that, to the first approximation, 
the field configuration in FR I jet over 10 kpc scales is a mix- 
tu re of toroidal and longit udinal components. 

iGizani & Leahy I d2004l) 's ROSATX-ray observations of the 
intracluster gas in the Here A cluster have revealed extended 
X-ray emission coming from a compact source at the center of 
cluster, out to a radius of 2.2 Mpc. The azimuthally-averaged 
X-ray surface bri ghtness profile is well fitted by a modified 
King model dKing I Tl972; Cavaliere & Fusco-Femiano 
119761) . with core radius R c = 121 ± 10 kpc and (3 = 
0.74 ± 0.03. Their X-ray analysis supports the conclusion 
that the inner jets are confined thermally by the cluster core 
gas: thermal pressure in the cluster core is almost 10 times 



larger than the equipartition pressure in the eastern jet with the 
standard assumptions (including no protons). The radio lobes 
are largely positioned beyond the X-ray core radius, allowing 
for projection, so they are expanding essentially into a power- 
law atmosphere with density falling as Br 313 ~ R~ 2 22 , quite 
close to t he R 2 profile needed to give the lobes a self-similar 
structure dFalle lll991l) . 

These observational evidence motivate us to apply our dy- 
namical magnetic tower jet model to Here A. To our knowl- 
edge, no MHD model has been applied to Here A before 
though there are several hydrodynamic models as we have 
already referred. Radio observations suggest that the mag- 
netic fields need to be taken into account to discuss about 
jet/lobe dynamics in Here A. We believe that a key clue to 
understanding the transition between the jet and the lobe may 
be whether or not the inner jets become over-pressured dur- 
ing their propagation in the X-ray thermal cluster gas. How- 
ever, the jet internal pressure may be magnetically, rather than 
thermally, dominat ed. This is in contrast to previous models 
dMeier et al. 1 199 lh . Of particular interest here are the jet/lobe 
transition, wiggling, and the magnetic field distribution asso- 
ciated with the eastern jet of Here A, which are not yet fully 
discussed among the past hydrodynamic models of Here A. 
We do not discuss the formation of rings in the western jet in 
this paper. 

This paper is organized as follows. In §2, we outline our 
numerical methods. In §3, we describe our numerical results. 
Discussions and conclusions are given in §4 and §5. 

2. NUMERICAL METHODS AND MODEL ASSUMPTIONS 

We solve the nonlinear system of time-dependent ideal 
MHD equations numerically in a 3-D Cartesian coordinate 
system (x, y, z). The basic numerical treatments (including 
the MHD numerical scheme) is essentially the same as that in 
Papers I - III. We assume an initial hydrostatic equilibrium in 
the gravitationally stratified medium, adopting an isothermal 
King model. The magnetic flux and the mass are steadily in- 
jected in a central small volume during a certain time period. 
Since the injected magnetic fields are not force-free, they will 
evolve as a "magnetic tower" and interact with the ambient 
medium. The dimensionless system of MHD equati ons is in- 
tegrat ed in time by using the TVD upwind scheme dLi & Li I 
120031) . Computations were performed on the parallel Linux 
clusters at the Los Alamos National Laboratory. 

2.1. Numerical Set Up 

We normalize physical quantities with the unit length scale 
Rq, the unit density po, and the sound speed C s o as the typical 
speed in the system. Other quantities are derived from their 
combinations, e.g., the typical time to is Rq/C s q, etc. In this 
paper, we use some observed quantiti es associated with the 
source Here A (IGizani & Leahy |[2004l hereafter GL04) in or- 
der to determine our normalization units. A central electron 
density n c = 1. Ox 10 -2 cm~ 3 is adoptedin GL04, suggesting 
quite a dense cluster, as such densities in clus ters are typically 
of order 10~ 3 cm -3 dJones & Forman IT9 84). This provides a 
unit density po = 1.7 x 10~ 26 g cm~ 3 . A single-temperature 
fit with kT = 4.25 keV is assumed for the /3-model cluster, 
giving a unit temperature To = 4.9 x 10 7 K, a unit pressure 
Po = 2.3 x 10~ 10 dyn cm -2 , and a unit velocity (sound speed) 
C s o(= 7Po/po) = 1-2 x 10 s cm s _1 . Here, we choose a unit 
length Ro = 30 kpc, and therefore, a unit time to = 2.4 x 10 7 
yr. The unit magnetic field B is (47rp C ,2 ) 1/2 = 53.4 pG. 
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In the King model we use here, we take the cluster core ra- 
dius i? c to be 120 kpc, corresponding to the normalized core 
radius 4.0, and take the exponential slope k to be 1.1 (3/3/2 
with (3 = 0.74) by adopting the observations in GL04. 

We use po and po as the initial quantities at the origin 
(x, y, z) = (0, 0, 0) and thus, the normalized p and p are set 
to unity at the origin in these simulations. The total compu- 
tational domain is taken to be \y\, \z\ < 15. The number 
of grid points in the simulations reported here is 360 3 , where 
the grid points are assigned uniformly in the x, y, and z direc- 
tions. The simulation domain is from —450 to 450 kpc, with 
Ax = Ay — Az ~ 0.083, which correspond to ~ 2.5 kpc. 

The injections of magnetic flux, mass, and the associated 
energies are the same as those described in Paper I. The ratio 
between the toroidal to poloidal fluxes of the injected fields 
is characterized by a parameter a — 25. The magnetic field 
injection rate is described by jj, and is set to be 7^ = 1. The 
mass is injected at a rate of j p = 0.1 over a central volume 
with a characteristic radius r p = 0.5. Magnetic fluxes and 
mass are continuously injected for = 3.0, after which the 
injection is turned off. These parameters correspond to a mag- 
netic energy injection rate of ~ 2.7 x 10 46 ergs s _1 and an 
injection time ~ 72 Myrs. The energy injection rate used here 
is a similar range to that of other hydro dynamic simulation s 
with the jet power 2 x 10 45 - 46 ergs s" 1 dSaxton et al. 1120021) . 
We use the outflow boundary conditions at all outer bound- 
aries. Note that for most of the simulation duration, the waves 
and magnetic fields stay within the simulation box, and all 
magnetic fields are self-sustained by their internal currents. 

2.2. Formulation of Synthetic Observation Images 

Based on our simulations, we make synthetic observation 
images to compare wi th the real observations of Here A 
dGizani & Leahy II2003L hereafter GL03). The synthetic ra- 
dio telescope observations presented here were produced us- 
ing a modified f orm of the syn t hetic observation technique 
first described in Tregilli s et al. I (1200 lb . Here we provide an 
overview of the key points relevant to the present investiga- 
tion; a complet e description of the technique is provided in 
lTregillisld2002h . 

Our method utilizes the vector magnetic field structures 
evolved self-consistently within these MHD calculations and 
a relativistic-electron momentum distribution, f(p), chosen 
by fiat during postprocessing. The function f(p) is defined 
over eight logarithmically-spaced momentum bins. Each bin 
has an associated slope value, q, and within each bin the dis- 
tribution is a power law (/ ~ p~ q )- Slope values may vary 
between bins, however, making it possible to specify non- 
power law distributions, such as would be appropriate for 
cases where radiative aging is significant. The specified dis- 
tribution functions are free to vary from location to location 
within the simulated flows, as is expected to be the case in 
real radio galaxies. We note that although these distributions 
are not derived self-consistently wit hin the computed flows, a s 
has been done in previous works dTregillis et al. 11200 lLl2004T) . 
imposing it as a postprocessing step does provide a measure 
of flexibility in our investigations. 

In this paper we choose the sign convention such that a > 
0, and the flux density S v cx v~ a . For the present investiga- 
tions, which do not incorporate the effects of radiative aging 
and reacceleration, we applied a spatially-uniform momentum 
slope value q = 4.5 throughout the simulated source. This 
corresponds to a synchrotron spectral index a = 0.75, which 
is consistent with the observed range of spectral indices within 



individ ual radio galaxies. For example, iGizani & Leahy I 
(120031) mapped the spectral index distribution in Here A be- 
tween 1.3 and 4.8 GHz, finding 0.6 < af 8 3 < 2.0, with the 
flattest values within the western jet and the steepest values 
along the lobe edges. They also found af ;§ « 0.75 in portions 
of the eastern jet. Our choice, a — 0.75, is also consistent 
with typical hotspot spect ral index values ° i n the sample 
of 3CR sources studied bv lAlexander & Leahy I (Q~987). 

The local number density of energetic electrons varied with 
the local magnetic energy density (i.e., n e cx B 2 ). This ap- 
proach represents a partitioning relationship between the field 
and particl e energies, suc h as is posited by minimum-energy 
arguments (lMileyl |1980). However, the overall scaling was 
chosen such that the energy density of relativistic electrons 
was far less than the magnetic energy density (i.e., we assume 
our simulated object is out of equipartition, which is consis- 
tent with the lack of dynamical feedback from the electrons in 
these calculations). 

Given information about the local magnetic field and the 
local distribution of energetic electrons, we compute a syn- 
chrotron emissiv ity j v in every zon e of the computational 
grid. As given by I Jones et al. I d 19741) . the emissivity is 

3 S (v)=J a0 ^ff(Ps)p q s ^f-yVB ± . (1) 

The spectral index a is related to the local electron momentum 
index q via a = (q — 3) /2, vb ± = eB sin £1/ (2Trm e c), where 
ft projects the local field onto the sky, and j a o is an order - 
unity dimensionless constant, defined in iJones et al. I d!974l) . 
For a selected observing frequency, v, the distribution, f{p s ), 
and the index, q, are determined for each point on the numeri- 
cal grid by establishing the relevant electron momentum from 
the relation^ = [2v/(3i'b_ l )]^ 1 ^ 2 \ withp s in units m e c. 

Once the synchrotron emissivity is known for every numer- 
ical zone, we compute surface brightness maps for optically- 
thin emission via line-of-sight integrations. Similarly, we 
can obtain the Stokes Q and U parameters for the syn- 
chrotron emission, as well as the correction for Faraday rota- 
tion through the source, making detailed polarimetric studies 
possible. We write the resultant data in FITS format, and an- 
alyze it usin g conventiona l observational packages (MIRIAD 
and KARMA (IGooch II 1 9951) 1 . 

3. RESULTS 

Our magnetic tower model has been applied to the Here A. 
In particular, our goal is to reproduce large scale structures in 
the jets: the jet/lobe transition, the non-axisymmetric defor- 
mation, and the projected magnetic field distribution associ- 
ated with the jets and lobes, as observed in GL03. 

Figure[T]shows a snapshot of the gas pressure p in the (x—z) 
cross section at t = 8.0 (192 Myr). As we have examined in 
papers II and III, the global structure of the magnetic tower 
consists of a well-collimated "body" and radially extended 
"lobes" in the gravitationally stratified atmosphere. A transi- 
tion from the narrow jet to the fat lobe can occur at the cluster 
core radius. This interesting property can be confirmed in re- 
sent observations (GL03, GL04). Several key features in these 
magnetic tower jets (the MHD wave structures, the heating 
process at the tower front, the cylindrical radial force equilib- 
rium at the tower edge, and the dynamic collimation process) 
are similar to the results of paper II (see also Fig. 1 of paper 
II for details of time evolution). The helically-twisted mag- 
netic tower has a closed current system that includes a pair of 
current circuits. Each circuit contains a forward electric cur- 
rent path (the jet flow itself, with its toroidal magnetic field, 
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4 5 6 7 

FIG. 1. — Distribution of the gas pressure p (logarithmic scale) with the 
velocity field (V x , V z ) (arrows) in the x — z plane at 192 Myr. The length 
and time scales are converted into the corresponding dimensional values. The 

cavities (low density and gas pressure) are formed due to the magnetic tower FIG. 2. — Axial profiles of physical quantities along the z-axis at t = 3.0 

expansion as corresponding lobes. The jets are deformed into the wiggled (t = 72 Myrs). Density p, sound speed C s , and the axial velocity component 

structures inside the lobes. V z . The positions of the expanding shock and lobe fronts are shown. 



toward the lobe) and a return current (along some path back 
to the AGN core). The global picture of a current-carryingjet 
with a closed current system linking the magnetosph ere of the 
centra l engi ne and the hot spots was introduced by iBenford I 
(119781 H006) and applied to AGN double radio sources. As 
seen in Fig. Q] the high-pressure material of the ambient gas 
seems to follow the edge of the jet near the bottom of each 
lobe. However, this is not the case, but the local pressure 
enhancement occurs inside the lobes. In the jet and lobe sys- 
tem, the axial currents drive the current-driven instabilities; 
it eventually produces magnetic reconnection process even in 
our ideal MHD assumption (the magnetic field will dissipate 
numerically). As a result, some local heating may occur, but 
it never affects the dynamics of jet propagation dramatically. 

For a full 3-D visualization of the magnetic lines of force, 
the reader is referred to paper III. The basic behavior of the 
helically-twisted magnetic field in the present paper appears 
similar to that in paper III. The magnetic tower jet has a well- 
ordered helical field configuration, with a tightly-wound cen- 
tral helix going up along the central axis and a loosely -wound 
helix coming back at the outer edge of the magnetic tower. 
The outer edge of the magnetic tower may be identified as 
a tangential discontinuity without the normal field compo- 
nent. The interior of tower (lobe) is separated from the non- 
magnetized external gas via this discontinuity. At the tower 
edge, the outward-directed magnetic pressure gradient force 
is roughly balanced with the inwardly-directed thermal pres- 
sure gradient force. On the other hand, at the core part of 
jet body, a quasi-force free equilibrium is achieved. In the 
context of magnetically-controlled fusion systems, the helical 
field in the magnetic tower can be regarded as the reversed 
field pinch (RFP) profile. 

The jet axial current and the ambient gas pressure can to- 
gether determine the radius of the magnetically-dominated 
lobes (paper II). That is, the internal gas pressure plays a mi- 



nor role in the lobes (even if th e thermal pressure is g reater 
than the non-thermal pressure (Morg anti et al. 1 119881) as i s 
typically seen in FR I radio galaxies dCroston et al. Il2003l) . 
The Alfven speed becomes large (about three times the local 
sound speed), while the plasma /3(= 2p/B 2 ) becomes small 
(f3 < 0.1) inside the density cavities due to the expansion of 
magnetic fluxes. As we have already discussed in paper II, it 
is remarkably noted that the expanding tower launches a pre- 
ceding hydrodynamic shock in the ICM, which may be asso- 
ciated with AGN-driven shocks seen in recent X-ray observa- 
tions dForman et al. Il200l iNulsen et"al~ll2005l; iFabian et al. I 
l2Q06h . 

Figure [2] exhibits physical quantities along the z-axis at 
t = 3.0 (t = 72 Myrs). The magnetic tower-driven hy- 
drodynamic shock front can be seen around z ~ 6.4 in the 
profiles of p, C s , and V z . The sound speed jumps by 27 % 
[the temperature T(oc C s 2 ) increases by 61 %] in the post- 
shocked region due to the shock compression. The magnetic 
tower front, which is identified as the lobe front, is located at 
z ~ 5.2. There is another MHD wave front at z ~ 4.5 where 
an increase in C s (density/pressure also) is accompanied by a 
increase in magnetic pressure (not shown in Fig. |2j, indicating 
the reverse MHD slow mode. In later time, this can steepen 
into a MHD slow shock that cause a compression (heating) at 
the head of magnetic tower as seen in Fig. [T](see also Fig. 2 
and §3.1 of paper II for details). 

Figure[3]displays the positions of the shock and lobe fronts 
as a function of time along the z-axis (propagating towards 
the positive direction). The shock front has a constant prop- 
agation speed Vsh ~ 2.1 (2.52 x 10 s cm s _1 ) that implies a 
Mach number ~ 1.63. On the other hand, the magnetically- 
dominated lobe front expands subsonically. This lobe front 
also has a quasi-constant propagation speed Vlb ~ 1-5 
(1.8 x 10 s cm s _1 ) that implies a Mach number ~ 0.91 (C s 
increases up to ~ 1.64 in shocked ICM; see also Fig. [2]). 
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FIG. 3. — The temporal evolution of the shock and lobe front positions 
along the z-axis. During the time evolution, both fronts propagate with al- 
most constant speeds. The shock front expands supersonically with Mach 
number ~ 1.63 (relative to the unshocked ICM), while the lobe front pro- 
ceeds subsonically with Mach number ~ 0.9 (relative to the shocked ICM 
cocoon). 



The current-carrying magnetic tower jet, which possesses 
a highly-wound helical field, is subject to the current-driven 
instability (CDI). Although the destabilization criteria will be 
modified by the ambient gas and the RFP configuration of 
the tower, we find that the propagating magnetic tower jets 
can develop the non-axisymmetric CDI modes. As seen in 
Fig. |U both the internal elliptical (to = 2) mode like the 
"double helix" and the external kink (to = 1) mode grow 
to produce the wiggles at different locations (see also pa- 
per III for details). We see no disruption owing to shear- 
driven, Kelvin-Helmholtz modes and/or the pressure-driven 
interchange modes. 

Figure [5] shows the synthetic radio synchrotron intensity at 
5 GHz, along with the polarization magnetic field vectors. 
The wiggled structures can be observed in both the jet and 
lobe. Also, filamentary structures appear beyond R c . The 
projected magnetic field is aligned with the jet and along the 
lobe edge. The magnetic tower model produces polarization 
features that are qualitatively consistent with those seen in 
the Here A observations of GL03. Furthermore, our dynam- 
ical solution of magnetic tower's evolution may give a self- 
consistent explanation to a feasible magnetic configuration 
discussed in non-dynamical models having "spine" (B^ > B z 
inside the lobe) - "Sheath" {B z > B^) structu r e in FR I jet 
over 1 kpc scale as discussed in Laing] (ll981h : ILaing et al~l 
(2006). 

The image contrast in our synthetic synchrotron surface- 
brightness map is approximately 2000: 1 between the brightest 
pixels (beams) near the central core and the faintest regions 
inside the lobes. Outside the small bright core region, the 
image contrast is roughly 200: 1 . 

When the synthetic 5 GHz synchrotron map is scaled to 
an angular size comparable to that of Here A, the brightest 
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FIG. 4. — Distribution of the axial current density J z in the x—z plane at the 
same time as Fig. fj] Two types of growing current-driven (external/internal) 
modes are visible in the lobe and jet. 




FIG. 5. — Synthetic synchrotron surface brightness computed at 5 GHz 
(greyscale), overlaid with polarization magnetic field vectors. The field vec- 
tors are corrected for Faraday rotation, and represent the field orientations 
at the source. The polarization field is nicely aligned with the jets and lobe 
edges. 



pixel (beam) in the synthetic 5 GHz surface brightness map 
is 875 mJy/beam. However, away from the small, bright re- 
gion near the core, the jets generally exhibit flux densities on 
the order of 30-70 mJy/beam, the central sheath exhibits flux 
densities in the range 20-30 mJy/beam, and the darkest parts 
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of the lobes exhibit < 1.0 mJy/beam. 

Our simulated source exhibits a fractional linear polariza- 
tion, m, at 5 GHz ranging from a few percent in the most 
magnetically tangled portions of the lobes to m *=s 30 — 50% 
on the lobe edges. These values are consistent with those 
typically observed in radio galaxies dBridle & Perly I [l984t 
iMuxlow & Garringtonlll99~ll) as well as with the m values re- 
ported in GL03 for Here A at 4848 MHz. Fractional polariza- 
tion values exceed 60% (roughly twice that reported in GL03) 
in small, isolated pockets in the sheath around the central 
jet and near the central core, which may indicate insufficient 
computational resolution near the core for properly captur- 
ing cellular depolarization. In regions outside the magnetized 
lobes, where the ambient medium in our calculational do- 
main has been mildly perturbed by expansion of the magnetic 
tower, the degree of polarization approaches 72%. This is 
the expected result for a region with a uniform magnetic field 
and a uniform p ower-law distribution of n onthermal electrons 
with q = 4.5 dRvbicki & Lightman|[l~979l) . 

4. DISCUSSION 

4.1. Shock and Lobe Expansions 

Some features about the shock and lobe expansions in the 
simulation are discussed by comparing with observations. A 
Chandra X-ray image of the Here A cluster shows that it has 
cavitie s and a shock front associated with the powerful radio 
lobes (Nul sen et al. 1 12005b . Unusually, these cavities show 
no clear connection to the radio sources, while many cavities 
are associated with radi o sources in other cluster cases (e.g., 
iMcNamara et al. 1120001). indicating that they might be ghost 
c avities (e. g., McNamar a et al. 11200 ll) . 

iNulsen et al. I d2005l) 's observation reveals that a shock 
front surrounds the Here A radio source. It is elongated in 
the dir ection of the r adio lobes and appears to be a cocoon 
shock dScheuerlll974l) . Their fitting of a simple hydrody- 
namic model to the surface brightness profile gives a Mach 
number for the shock front of ~ 1.65 and its total energy 
~ 3 x 10 61 ergs (its mean mechanical power ~ 1.6 X 10 46 ergs 
s _1 ) is estimated. As we examined in Sj3] the magnetic tower 
drives a hydrodynamics shock front in the ambient ICM. De- 
rived Mach number M ~ 1.63 matches quite well to their 
observation and also our energy inputs are in qualitatively rea- 
sonable range. 

As the radio lobes of Here A lack anything resembling 
hotspots, it seems likely that they are expanding at around 
the sound speed, and so are confined by the thermal pres- 
sure of the ambient I CM, rather than by strong shocks 
(lGizani& Leahy 1 12004). To our knowledge, there may be 
no direct result revealing the lobe expansion speed of Here 
A radio lobes, but we discuss general properties o f expanding 
lobes in radio galaxy based on Kra ft et aLl (|2006)'s argument. 
It is believed that the radio lobes are greatly overpressurized 
relative to the ambient medium, indicating that sharp X-ray 
surface brightness discontinuities or large jumps in gas tem- 
perature at the lobe edges. A high Mach number (M ~ 8.5) 
shock is detected around the SW radio lobe of Centaurus A 
dKraft et al. Il2003l) . 

However, there is no evidence especially for FR I 
or low-power FR II radio galaxie s, including Her e A 
dGizani & Leahy 1 12001 and 3C 388 dKraft et al l 1 120061) . If 
thermal conduction in the ICM is efficient (an order of the 
Spitzer value), the n shocks may be near ly isothermal and thus 
difficult to detect (Fabia net al. Il2006h . However, Chandra 
observations suggest that the thermal conduction (and viscos- 



ity) o f the ICM is orders of magnitude below th e Spitzer value 
(e.g.. iVikhlinin et al. Il200ll: iKraft et al. Il2~003l Thus, some 
lobes including Here A may be inflating transonically and 
subsonically (i.e., M < 1). Based on hyd rodynamic argu- 
ments of a buoyant bubble, M ~ 0. 6 in M87 ( Churazo v et al.1 
|2002J) and M ~ 0.5-0.9 in 3C 388 dKraft et al. 120061) are es- 
timated. Our numerical result indicates that the magnetically- 
dominated lobe expands subsonically with M ~ 0.91, dis- 
cussed in ij3] This will also match the observationally prefer- 
able picture. 

4.2. Lobe Size 

We next dis cuss the radial size of the lobes n be as based 
on our model dNakamura et al. f2006l) . The poloidal magnetic 
fluxes and the poloidal current I z in the lobes remain well 
collimated around the central axis even after the system has 
evolved fully. This implies that the toroidal magnetic fields in 
the lobe region are distributed roughly as 



h/r. 



(2) 



When the lobe experiences sufficient expansion, we can ex- 
pect \B Z \ ~ \BA at the outer edge of the RFP configuration 
dFreidbergl 1 1 9821) . The magnetic and background pressures 
tend to balance each other at the tower edge. Thus, we can 
expect that 

Bl + B\ 

-~B\ 



J 4> 



(3) 



where p c is the external gas pressure at the tower edge. By 
combining eqs. (O and OJ, we have 



B 



I, 



Hobc 



which gives 



Hobo 



p, 



-1/2 



(4) 



(5) 



Figure [6] shows the distribution of p along a transverse line 
at z — 8.5 (255 kpc), when the lobe has a maximum radius 
(after 192 Myr, as in Figs. Q]and[4]i. We can see that the exter- 
nal gas pressure begins to decrease around x ~ 5 (150 kpc) 
toward the central (z) axis. The pressure cavity corresponds to 
the lobe of our magnetic tower. We can estimate the lobe ra- 
dius Hobe predicted by eq. © from physical quantities in sim- 
ulation. From Fig. we take J z ~ 3 (4.5 x 10~ 24 A/cm" 2 ) 
within a radius r ~ 0.75 (22.5 kpc), which makes I z ~ 1.68 
( 5 x 10 17 A). As suggested by Fig. [4] we estimate p c ~ 0.11 
(2.5 x 10 -11 dyn cm -2 ) for the external gas pressure. We 
therefore obtain n bc ~ 5.2 or roughly 156 kpc. Thus, we 
find that the depression radius in Fig. [6] is qualitatively con- 
sistent with n be and is c omparable with the phy sical scale 
seen in radio observations dGizani & Leahy 1 12003). Based on 
this numerical result, we estimate that one simulated lobe of 
Herc-A may contain a total energy ~ 1.8 x 10 61 ergs, out of 
which ~ 1.0 x 10 61 ergs are in magnetic energy, ~ 1.4x 10 60 
ergs are in kinetic energy, and ~ 6.6 x 10 60 ergs are in ther- 
mal energy, respectively. In order to explain the lobe radius 
by this magnetically dominated lobe system, an axial current 
~ 5 x 10 A, which is flowing along the jet central axis, is 
needed. 

5. CONCLUSIONS 

By performing 3D MHD simulations, we have investigated 
the nonlinear dynamics of magnetic tower jets in a galaxy 
cluster environment. In this fourth paper of the series, we 
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FIG. 6. — Transverse profile in the re-direction of the gas pressure p on 
z = 8.5 (255 kpc). The lobe radius riobc calculated via eq. {5) is also 
shown. The pressure cavity radius and ri t, c qualitatively agree with each 
other. 



have applied our numerical modeling of the magnetic tower 
to a specific source, Hercules A. 

To our knowledge, this is the first trial of discussing Her- 
cules A jet by using a magnetohydrodynamic model. Pro- 
jected magnetic field distributions associated with the jets and 
lobes of Hercules A have been revealed in recent radio obser- 
vation and therefore, the dynamics including magnetic fields 
plays a role in understanding this source. In the present pa- 
per, we have investigated the jet-driven shock, jet/lobe tran- 
sition, wiggling, and projected magnetic field distribution of 
Hercules A and, in general, these feasutres can be explained 
by our magnetic tower model. 

Our conclusions are summarized as follows: 



2. The expanding magnetic tower jet produces the pre- 
ceding hydrodynamic shock wave in the ambient ICM. 
This can be interpreted as AGN-driven shock in X- 
ray observatio n and the derived shock Mach number 
is identical to lNulsen et al. I d2005)'s observation result. 
The magnetic tower lobes expand subsonically and thus 
no hot spots are produced at the end parts of lobes. This 
may be a general un derstanding of FR I or low power 
FR II radio galaxies (iKraft et al. 1120061) . 

3. The size of the magnetic tower lobe can be accurately 
estimated from the jet axial current and the external gas 
pressure at the lobe edge. The estimated lobe size de- 
rived by our numerical simulation is qu alitatively com- 
parabl e with the physical scale seen in lGizani & Leahvl 
d2003l In our model, if the background gas has a steep 
slope (fi) in the King profile, the magnetic jet can po- 
tentially lead to a huge magnetic tower lobe. 

4. Magnetic tower jets, which have a reversed magnetic 
pinch profile, are subject to current-driven instabili- 
ties. Specifically, the apparent distortions by the non- 
axisymmetric kink mode are visible at the jet body in- 
side the lobe rather than the central core. By expanding 
the lobe the edge of axial currents could be free against 
the exciting external kink mode; this situation may be 
suppressed inside the cluster core radius. 

5. The synthetic polarimetry of our magnetic-tower jets is 
consistent with the gross polarization features observed 
in the jets of Hercules A. The magnetic tower model 
produces projected B-vector distributions that are sim- 
ilar to those observed in the global extragalactic jet-lobe 
system, especially along the jets and the lobe edges. 

6. Magnetic tower jets may produce a self-consistent pic- 
ture of the field configuration discussed in several non- 
dynamical model of FR I jets. It consists of toroidally 
dominated spine inside the lobe and poloidally domi- 
nated sheath at the lobe edge (a closed poloidal mag- 
netic flux with a toroidal magnetic component, rather 
than one side directed helical configuration). 



1. A transition from the narrow jet to the fat lobe can 
occur at the cluster core radius. This interesting 
feature can also be confirmed in recent observations 
dGizani&Leahvl I2003L |2004|) . In our models, the 
tightly collimated helical field configuration of a mag- 
netic tower jet expands abruptly beyond the core radius 
due to the decreasing external gas pressure. Thus, the 
pressure profile of the surrounding ICM plays an impor- 
tant role in determining the morphology of a magnetic 
tower jet. 
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